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Abstract. We present a matched- filter based algorithm for transit detection 
and its application to simulated COROT light curves. This algorithm stems 
from the work by Borde, Rouan & Leger (2003, A&A 405, 1137). We describe 
the different steps we intend to take to discriminate between planets and stellar 
companions using the three photometric bands provided by COROT. These steps 
include the search for secondary transits, the search for ellipsoidal variability, 
and the study of transit chromaticity. We also discuss the performance of this 
approach in the context of blind tests organized inside the COROT exoplanet 
consortium. 



1. Introduction 



COROT (BagHn][2003), scheduled for launch by the end of 2006, will be the 



first space mission with the cap ability to detect ex trasolar planets with sizes 
down to a couple of Earth radii ( Rouan et al. |[2000l ). It will observe a total of 
60,000 stars in 5 runs of 150 days each. In a previous paper, we worked out a 



theor etical estimate of the planet detection capability of COROT (jBorde et al. 
2003). In this paper, we first describe our transit detection algorithm, dubbed 



Transdet, then we report its performance on simulated light curves (LCs), and 
finally, we discuss the tools we are developing to discriminate between genuine 
exoplanets and astrophysical false positives (binary stars). 
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2. Key ideas about Transdet 

At low signal-to-noise (SNR), a transit signal can be approximated by a rectan- 
gular signal oscillating between a high level H (outside the transit) and a low 
level L (during the transit). The amplitude of this signal depends on the ratio 
of the planet surface to that of the star: H — L = eH, where e = {R p / R*) 2 . The 
time spent at level L, or transit duration r, is much shorter than the transit 
period T. We note to the date of the beginning of the first transit falling inside 
the observation window (150 days for COROT). 

In order to detect the transit signal in a noisy LC, we adopt the matched- 
filter approach that consists in correlating the LC with rectangular signals of 
the type described above. We call test-signals these rectangular signals. Each of 
them is characterized by a parameter triplet (T, f , to)- The correlation reaches a 
peak when (T, f , to) = (T, t, to). Practically LCs must be first high-pass filtered 
to remove irrelevant low-frequency variations (stellar fluctuations) that would 
bias the correlation, and we take the negative of the LCs so that transits would 
produce positive correlation peaks. 

For, this method to work, the three-dimensional space (T, f,to) must be 
explored with a sufficient resolution in the given range accessible to COROT 
(1 < T < 75 days, 1 < f < 10 hrs, < to < T). The choice of this resolu- 
tion results from a trade-off between computing time and detection efficiency. 
Typically, we use steps (5T, 5f, <5to) small enough so that the maximum corre- 
lation signal would be no less than 75 % of its theoretical maximum value. This 
requirement translates into the computation of 14 x 10 6 correlations per LC, a 
task completed in 90 s with a fully IDL-coded algorithm on a Pentium M at 
1.6 Ghz. However, preliminary tests indicate that an optimized IDL-C hybrid 
code could cut down the computing time by a factor 7-8, which would make 
possible a higher detection efficiency. 



3. Detection threshold 

3.1. Theoretical detection threshold 

Practically, we use as our primary detection metric not the correlation itself, 
but the SNR on the correlation defined as 

SNR% = (fcr) 1 / 2 eH/a, (1) 

where a is the standard deviation of the noise affecting the high-pass filtered 
LC, and k is the number of transits in the observation window. 

A theoretical detection threshold can be simply set assuming: (1) a total 
number of L Cs, (2) a false detection rate, and (3) Gaussian white noise (see 



Borde II2003I . p. 45). Requiring less than 1 false detection, we obtain SNRi > 6.5 



for 1,000 LCs, and SNRi > 7.0 for 60,000 LCs. 
3.2. Empirical detection threshold 

In order to train the co-investigators of the COROT exoplanet consortium and 
compare the performances of existing transit detection algorithms, a first tran- 
sit detection blind test (BT1) was organized inside the consortium in 2004 



Transdet: a transit detection algorithm 



3 



( Moutou et al. 2005). BT1 consisted in the analysis of 999 LCs with no in- 
formation on the target stars. 

Figure Q] (left) shows the distribution of SNR\ for BT1. A priori, one would 
expect a bimodal distribution, where LCs without transits would appear at low 
SNR\, whereas LCs with transits would appear at high SNR\, and a dividing 
line around 6.5. If the obtained distribution is indeed bimodal, the plot suggests 
a dividing line, or detection threshold, of either « 15 (175 selected LCs) or w 30 
(28 selected LCs), largely above the theoretical estimate. This discre pancy stems 



from t he inaccuracy of the white noise hypothesis as pointed out by I Pont et al. 
(|2006h . 



BTl: SNRl histogram (999 LCs) 



BT1: SNR2 histogram (999 LCs) 





Figure 1.: Left: SNRi histogram. The dashed and dotted lines mark thresholds 
at 15 and 30, respectively. Right: SNR2 histogram. The dashed and dotted 
lines mark thresholds at 4 and 5, respectively. 

We find that a secondary detection metric computed on the folded LC 
(FLC) helps the selection process. We define this new metric as 



SNR 2 = (H-L)/a f , 



(2) 



where H and L are measured on the FLC, and ay is the standard deviation 
of the noise affecting the high-pass filtered FLC. Figure [T] (right) shows the 
distribution of SNR2 for BTl. This new plot suggests a threshold of either ~ 4 
(42 selected LCs) or « 5 (16 selected LCs). It turns out that SNR 2 > 4.0 is the 
right criterion as the 12 genuine planets correctly detected by Transdet in this 
exercise are among the 42 selected LCs (see also Fig. [2]). 



3.3. Detection capability 

Figure [3] is an attempt to estimate the detectio n capability of Trans det from 
the results of BTl in much the same way as in iMoutou et aTl (|2005l . Fig. 7). 
The solid line is a line at constant SNR± adjusted so that all detected planets 
lie above it, and all undetected planets lie under it. From this, we conclude that 
the detection limit would be about 300 ppm at r/T = 0.04, which could be a 2 
Earth-radii planet around a G2V star with a 3-day period. 
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BTl: SNR2 vs. SNR1 (999 LCs) 
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Figure 2.: SNR2 vs. SNRi diagram. Crosses and squares are overplotted on 
genuine planets and binary stars, respectively. A thresholding in SNR2 selects 
all planets much more efficiently than a thresholding in SNR\. 



BTl: Transdet detection capability 
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Figure 3.: Comparison between actual and measured values for e and t/T for 
BTl. A line at constant SNR\ marks the bottom of the region where transits 
were detected. 



4. False positive discrimination 

The vast majority of transit signals ar e not caused by planets, but rather by 
stellar companions in binary stars (e.g.. [Brown |[20(KS): grazing eclipsing binary 
(GEB) or background eclipsing binary (BEB) blended with foreground stars 
may produce transits with amplitudes consistent with planetary companions. 
Radial velocimetry provides the definitive way to discriminate planets from bi- 
naries. However, it would be much too time-consuming to observe every single 
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candidate, thus it is mandatory for the success of transits searches to sort out 
as many candidates as possible using the photometry alone. A second transit 
detection blind test (BT2) was organized inside the COROT consortium in 2006 
to work toward this goal. BT2 participants were given 237 LCs measured in the 
three COROT photometric channels, as well as the spectral types of the target 
stars. 

4.1. Anomalies on folded light curve 

A first clear indicator of the stellar nature of the companion is the presence of a 
secondary eclipse in the LC. Therefore, once the primary eclipse has been found, 
we scan the FLC for shallower transit-like events. 

A second indicator is the presence of ellipsoidal variability: tidal forces 
exerted by a massive companion at short orbital distance distort the primary's 
shape into an ellipsoid. Because the ellipsoid rotates with the orbital motion, 
the apparent stellar d isk varies and the LC's ba seline is modulated at twice the 
orbital frequency (see lSirko Sz Paczvhski 1 12003! ) . We also scan the FLC for this 
modulation. 

4.2. Transit parameter check 

COROT's entry catalog contains photometrically-determined spectral types for 
target stars, as well as estimates for stellar masses and radii. With the target 
radius and the measured transit depth, it is straightforward to check that the 
companion's radius lies in the planet range, say R p < 1.5 Rj. 

For short orbital periods where the assumption of a circular orbit is reason- 
able, the transit duration can be checked against the maximum transit duration 
computed for an equatorial transit, i.e., r max = 1.8 -R^T/M*) 1 ' 3 , where r max is 
in hours, T in days, and R* and M* in solar units. A measured duration greater 
than r max would either point to a problem with the transit parameters (e.g., an 
incorrect period) or with the parameters of the star in front of which the transit 
is occurring, thus betraying a BEB with a spectral type different than that of 
the target. 

Finally, the transit shape can be an indicator of GBs for partial occultations 
produce V-shape transits. To quantify this, we fit a trapeze to the transit in order 
to measure the duration Tp of the flat part of the transit (between ingress and 
egress) as well as the total transit duration r, and we compute VS = 1 — rp/r. 
Currently, our V-shape criterion is arbitrarily defined as VS > 0.8. 

4.3. Transit chromaticity 

COROT will have the unique capability to measure transits in three spatially- 
separated color channels thanks to a prism objective. Only BEBs perfectly 
aligned with target stars will produce transits in the three channels. Therefore, 
transit chromaticity will be a powerful indicator of a false positive. For the sake 
of simplicity, COROT's color channels will be referred to as B (blue), G (green), 
and R (red), but note that they do not correspond to any existing photometric 
systems. We define now three chromaticity indicators, 

K + 4) 1/2 K + 4) 1/2 (4 + 4) 1/2 ' [) 
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where eb, £g> and £r are the relative transits depths in the color channels, and 
as, &g, and an, the corresponding standard deviations of the noise. Figured] 
shows chromaticity diagrams for a selection of BT2 LCs. In this case, we applied 
a 3 cr threshold to classify LCs as chromatic. 
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Figure 4.: Chromaticity diagrams for a selection of BT2 LCs. Dotted boxes mark 
the 3<t limit. 



5. Conclusion 

We have developed Transdet, a matched-filter based algorithm for transit detec- 
tion, as well as tools for false positive discrimination using the LCs alone. All of 
these were developed and tested using blind tests conducted within the COROT 
consortium. However, fine tuning of the thresholds for the various tests will be 
done by studying the statistics of real COROT data. 
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